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ABSTRACT 

We use the Millennium Simulation series to study how the dynamical state of dark 
matter halos affects the relation between mass and concentration. We find that a large 
fraction of massive systems are identified when they are substantially out of equilib- 
rium and in a particular phase of their dynamical evolution: the more massive the 
halo, the more likely it is found at a transient stage of high concentration. This state 
reflects the recent assembly of massive halos and corresponds to the first pericentric 
passage of recently-accreted material when, before virialization, the kinetic and poten- 
tial energies reach maximum and minimum values, respectively. This result explains 
the puzzling upturn in the mass-concentration relation reported in recent work for 
massive halos; indeed, the upturn disappears when only dynamically-relaxed systems 
are considered in the analysis. Our results warn against applying simple equilibrium 
models to describe the structure of rare, massive galaxy clusters and urges caution 
when extrapolating scaling laws calibrated on lower-mass systems, where such de- 
viations from equilibrium are less common. The evolving dynamical state of galaxy 
clusters ought to be carefully taken into account if cluster studies are to provide precise 
cosmological constraints. 
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1 INTRODUCTION 

Galaxy clusters are powerful cosmological probes. Because 
of their large masses, they can be detected at cosmological 
distances despite being rare and studied at various wave- 
lengths, from X-rays to optical to millimeter wavelengths. 
Their rarity is actually a strength when it comes to cos- 
mology; massive clusters trace the tail of high-mass objects 
able to collapse under their own gravity and therefore their 
numbers, as a function of mass and redshift, are exponen- 
tially sensitive to the normalization of the power spectrum of 
density fluctuations and to the cosmological par ameters that 
gover n the universal expansion history (see, e.g.. lAllen et al.l 
I2OIII . for a recent review and a complete list of references). 

Cluster counts are thus widely recognized as a pre- 
mier tool able to provide cosmological constraints com- 
plementary to those based on analysis of the cosmic mi- 
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crowave background, supernova luminosity dist ances, galaxy 
clustering, and gravitat i onal lensing (see, e.g.. iMantz et all 
\ 'Cunha ct al.l 2009l : iHenrv et all2009l : IRozo et al.ll2009l : 
)0i). Thei 



Vikhlinin ct al. 20091). Their potential is more directly real- 
ized, however, when observables can be turned into effective 
measures of cluster mass; indeed, the most robust and dis- 
criminating theoretical predictions concern the abundance 
and redshift evolution of clusters of different mass. 

Mass, however, is not directly observable, which means 
in practice that observables such as cluster richness, veloc- 
ity dispersion, X-ray luminosity /temperature, or Sunyaev- 
Zel'dovich (SZ) decrement must be related to mass through 
scaling laws whose shape and scatter need careful calibra- 
tion. Although virial equilibrium and self-similarity suggest 
ower-law scalings between mass and observational proxies 
Kaiseij[l98^ :lN avarro et al.ll 19951 ). it has long been appreci- 
ated that baryon physics breaks self-similarity and imposes 
scalings with different slopes and re dshift dependence (e.g., 
lKaiseilll99ll : lEvrard fc Henrvlll99ll ') . 
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The complexity of the baryon physics involved pre- 
cludes robust theoretical predictions and therefore the mass- 
observable relations are usually parameterized as power laws 
of adjustable slope, redshift-dependent normalization, and 
Gaussian scatter that are empirically calibrated using a 
(usually small) set of well-studied clusters. Besides simplic- 
ity, the underlying rationale of this procedure is that, how- 
ever intricate, the mass-dependence of the physical processes 
responsible for a given mass-observable relation is monotonic 
and that deviations of individual clusters from the mean 
trends are driven by stochastic effects. 

The validity of these assumptions, however, is not as- 
sured and recent work has highlighted their limitations. One 
example is the mismatch between the average SZ signal at 
given optical richness measured by the Planck satellite and 
expectations based on the X-ray properti es of clusters with 
weak lensing-calibrated masses (^Planck Collaboration et al.l 
This might either indicate a dichotomy in the gas 
content of clusters of different mass or perhaps arises from 
the large, correlated scatter betwee n observables sugge sted 
by recent simulation work (see, e.g.. lAngulo et al]|2012l ). 

A second example is provided by recent reports of a non- 
monotonic relation betwee n the mass and concentration of 
mass ive dark matter halos (|Klvpin et al.|[201ll : IPrada et al.l 
[20l3). These authors report an "upturn" in the median con- 
centration on mass scales where the rms fluctuations of the 
linear density field, cr(M, z) ^ 1/2. This corresponds to rare 
objects with virial masses equivalent to a few times 10^* Mq 
at 2: = 0, shifting to lower masses at higher redshift. 

The origin of the upturn is unclear. Since halo con- 
centration reflect s the mean density of the universe at the 
time of assembly jNavarro et al.|[l996l,ll997l:[Eke et al.ll200ll : 
iBuUock et al.ll200ll : IWechsler et aLll2002l : IZhao et al.ll2003l) . 
the mass-concentration relation is expected to be a mono- 
tonic one where concentration decreases with increasing halo 
mass. An upturn in concentration may have non-negligible 
consequences on cluster studies: concentration is a mea- 
sure of the characteristic density of a cluster, and it could 
therefore affect directly mass proxies that are sensitive to 
density, such as X-ray luminosity. A non-monotonic mass- 
concentration relation might thus lead to selection biases or 
"breaks" in the relation between X-ray properties and mass 
that would be poorly captured by the assumed power-law 
scalings. 

We use here data from the Millennium Simulation se- 
ries to revisit the mass-concentration relation of rare, mas- 
sive halos. In particular, we examine the physical origin of 
the reported "upturn" in concentration and the nature of 
the scatter about the mean relation. We describe briefly the 
numerical simulations in Sec. [2l present our results in Sec. [3l 
and summarize our main conclusions in Sec. [D 



2 NUMERICAL SIMULATIONS 

2.1 The Millennium Simulation series 

We use in our analysis halos identified in the three "Mil- 
lennium Simulations", which we refer to as MS-I, MS-II, 
and MS-XXL, or collectively as MS. All MS runs adopted 
the same cosmology (a WMAP-1 normalized fiat ACDM 
model) and the same sequence of outputs times in order 



to facilitate comparisons between the runs. The values of 
the cosmological parameters are as follows: Sim ~ 0.25, 
r^A = 1 - r2m = 0.75, h = 0.73, n = 1 and ag = 0.9. 
Here fii is the present-day contribution of component i to 
the universal matter energy density in units of the critical 
density for closure; ag is the rms mass fiuctuation in spheres 
of 8h~^ Mpc radius linearly extrapolated to 2; = 0; n is the 
spectral index of primordial density fluctuations, and h is the 
Hubble parameter defined so that Ho = H{z = 0) = 100 h 
km/s/Mpc. 

The MS-I used 2160^ particles of mass iTir, = 8.61 x 
10** h~^MQ to follow the evolution of the dark mat- 
ter component in a periodic box 500 /i~^Mpc on a side 
(|Springel et al.ll2005l ). MS-II used the same number of par- 
ticles, but a box size 5 time s smaller, Lbox = 100 /i~^Mpc 
(|Bovlan-Kolchin et al.l l2009l l: each MS-I particle is thus 
125 X heavier than in MS-II. The MS-XXL is substantially 
larger than MS-I, both in box size (Lbox = 3/i~^Gpc) and 
total particle number (6720"*), and thus provides the best 
statistics for the rarest and most massive objects at any 
redshift (|Angulo et al.ll2012l ). albeit with poorer mass res- 
olution. The particle mass in the XXL is nip = 6.17 x 
10^/i~^Mq. The Plummer-equivalent gravitational soften- 
ing lengths used in these runs are tp = 5/i~^kpc (MS-I), 
1/i-^kpc (MS-II), and 10/i^^kpc (MS-XXL). 



2.2 Analysis 

Halos in the MS a re identified usin g a friends-of-friends 
(FOF) group finder (|Davis et al.lll985l ') with a linking length 
equal to 20% of the mean interparticle separation. Sub- 
structure in each FO F group is identified using SUBFIND 
(jSpringel et al.|[200lh . We consider in our analysis only the 
main halo of each FOF grouping, and define its virial mass, 
M200, as the mass contained within a sphere centered at the 
potential minimum that encloses a mean overdensity of 200 
times the critical density for closure, pcrit = SH{z)^ /8nG. 
This implicitly defines the virial radius of a halo, r^oo, and 
its virial velocity, V200 = y/ GM200 /'■200 • 

We adopt two different quantities to estimate halo con- 
centrations, c. One uses the ratio of the peak circula r speed , 
Vinax, to the virial velocity, V200, as in iPrada et all (|201lf ). 
Concentrations are derived from Vmax/V2oo assuming that 
the mass p rofile follows the NFW formula jNavarro et al.l 
1 19961 , [1993). We use all particles within r2oo to estimate 
c; i.e., Knax is computed without subtracting substructure. 
This is important in order to ensure convergence on the over- 
lapping mass scales of different MS runs because substruc- 
ture, which may affect the estimates, is particularly sensitive 
to numerical resolutiorQ. 

Some advantages of this estimator include the fact that 
it is simple to compute, fairly stable numerically (Knax is 
well defined and robustly estimated even at modest numeri- 
cal resolution) , and that it does not rely on profile fitting and 
residual-minimization techniques. The main disadvantage is 
that, because it assumes an NFW profile, it is unduly sen- 
sitive when Vmax approaches V200: for example, when con- 



^ The slight co ncentration mismat ch betwen MS-I and MS-II ha- 
los reported bv lPrada et al.l 1I2OIII) is in all likelihood duo to this 
effect (see Fig. [TJ. 
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Figure 1. Mass-concentration relation for all halos resolved with more than 5000 particles in the Millennium Simulations. The top 
panels show, as a function of virial mass, concentrations estimated via Vmax/V200i the ratio between the maximum circular velocity and 
the virial velocity of a halo. Bottom panels show concentrations estimated by the ratio between the half-mass radius and virial radius, 
Tioo/r^- For each set of panels, the corresponding NFW concentration values are indicated by the right-hand side tickmarks. Panels on 
the left correspond to halos identified at 2;^ = 0, those on the right at 2;^ = 3. The three groups of points in each panel correspond, from 
left to right, to MS-II, MS-I and MS-XXL halos, respectively. All halos are shown in grey, relaxed halos in orange. Connected symbols 
(see legends) indicate the median concentration computed in bins of virial mass. Open and filled symbols correspond to "relaxed" or all 
halos, respectively. Note the excellent agreement in the median concentration of halos of similar mass identified in different simulations. 
The bottom tickmarks indicate halo virial mass; those at the top the dimensionless mass measure, u = <5crit/'''(Af, z). Note that, because 
of the fixed comoving box size of each run, halos identified at 2 = 3 sample "rarer" peaks in the mass spectrum (i.e., higher values of u) 
than those at 2 = 0. 



centrations double from c = 2.5 to 5 V;nax/V200 changes 
only from 1.001 and 1.062, respectively. Even tiny changes 
in Knax may thus lead to large variations in the derived 
concentration. Further, a large value of l^max/1^200 does not 
necessarily imply a high central concentration of matter. For 
example, a massive substructure may elevate V^nax but, if lo- 
cated at the outskirts of the halo, it might actually decrease 
the average central density of a halo. 

Because of this, we adopt an additional concentration 
estimator: the virial-to-half mass radius ratio, r2oo/''h, pro- 
vides a more reliable tracer of the concentration of dark 
matter. As for the velocity ratio, we use all particles within 
'"200 to compute rh and express the ratio as a concentration, 
c, assuming an NFW profile. 

Our analysis considers all halos resolved with at least 



5000 particles (i.e., A^20o S5 5000) at three different red- 
shifts, 0, 1, and 3. This corresponds to a minimum mass 
of 3.44 X 10^" R-^Mq for MS-II, 4.30 x 10^^ R-^Mq for MS- 
I, and 3.09 x IQ^^h'^MQ for MS-XXL. There is a fairly large 
overlap in the mass scales covered by the different simula- 
tions, which allows us to check the sensitivity of our results 
to numerical resolution. 

As discussed bv lNeto et all (|2007l ). halo concentrations 
can be affected by transient departures from equilibrium, 
resulting in subtle biases in their mean mass depen- 
dence and scatter. We therefore identify a subsample of 
"relaxed" halos, as those that satisfy the following three cri- 
teria: (i) /sub < 0.1, (ii) doff < 0.07 and (iii) 2T l\U\ < 1.3. 
Here /sub = Msub/M2oo is the mass fraction contributed by 
substructure; doff ~ |rp — rcM|/f2oo is the offset between the 



4 Ludlow et al. 



position of the potential minimum and the halo barycenter, 
expressed in units of the virial radius; and 2T /\U\ is the 
virial ratio of kinetic to potential energies. 

A full discu ssion of these criteria is provided by 
iNeto et all (|2007l '). who point out the need for multiple ob- 
jective criteria to reliably assess the equilibrium state of a 
halo. This is because all of these measures fluctuate dur- 
ing the virialization process and may therefore briefly fail 
to identify out-of-equilibrium systems. Ifowever, the fluctu- 
ations of the three different measures are not synchronized 
and, therefore, it is unlikely that all three will fail at the 
same time. For example, a halo may have a large abundance 
of substructure but a small doff if its subhalos happen to 
be distributed approximately isotropically about the halo 
center. Conversely, /sub may be small and doff large when 
deviations from equilibrium are driven by a relatively small 
fraction of the mass, such as in a succession of minor merg- 
ers. Finally, doff and /aub make no use of dynamical informa- 
tion and therefore they often fail at weeding out particular 
phases of ongoing mergers, for example, when the merging 
partner is at pericentric passage. This configuration mini- 
mizes doff (most of the mass is near the center) and also 
/sub (substructures are more difficult to identify when they 
are close to the center) , but their transient nature would lead 
to large values of the virial ratio, 2T/\U\. The total number 
of halos considered in each simulation, as well as the to- 
tal number of relaxed halos, are listed for each redshift in 
Table [T] 



3 RESULTS 

3.1 Mass-concentration relation 

Figure [T] shows, as a function of virial mass, the concentra- 
tion of halos in our sample, as measured by Vmax/V2oo (top 
panels) or r2oo/^h (bottom panels). Tickmarks on the right 
express these quantities in terms of NFW-based concentra- 
tion values. All individual halos are shown by grey dots; 
those that pass the three relaxation criteria are shown in 
orange. Panels on the left show halos identified at Zid = 0; 
those on the right correspond to Zid = 3. Median concen- 
trations for all (blue) and relaxed (red) halos are shown by 
connected symbols (see legend on figure for identification). 
The excellent agreement between the different runs on over- 
lapping mass scales indicates that our results are unaffected 
by numerical resolution effects. 

The results for all halos shown in Fig. [1] confirm the up- 
tu rn in the med i an co ncen tration of massive halos reported 
bv lKlvpin et all (|201ll ) and lPrada et all (|201ll '). As discussed 
by these authors, the upturn affects only rare, extremely 
massive systems, and is therefore more easily appreciated at 
earlier times in simulations of fixed comoving box size (i.e., 
Zid = 3 in Fig. [l]), when massive halos trace rarer peaks 
of the density fluctuation field. We quantify this using the 
dimensionless "peak height" mass parameter 



y{M,z) = 5,,it{z)/G{M,z), 



(1) 



defined as the ratio between the critical overdensity for col- 
lapse at redshift z and the linear rms fluctuation at z within 
spheres containing mass M. The larger the value of the 
rarer the halo and the more massive it is relative to the char- 
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Figure 2. An illustrative example of how concentrations vary 
during the virialization process that follows a major episode of ac- 
cretion. The upper panel shows, as a function of redshift, the evo- 
lution of the virial mass, M200, and of the half-m ass radius con- 
centr ation estimator, 7'20o/''hi for Aquarius halo F l lSpringel et al.l 
This system, which has a half-mass formation redshift of 
= 0.56, is in equilibrium at 2 = according to all relaxation 
criteria, whose evolution is plotted in the bottom two panels. 
Large symbols are used to indicate times when the halo would 
not be considered relaxed, according to at least one of the re- 
laxation criteria. Although the halo undergoes a major merger 
at 2 0.43, by 2 = more than a crossing time has elapsed, 
allowing the system to reach equilibrium. During virialization, 
however, the concentration estimator, r20o/i"h, fluctuates by al- 
most a factor of two. It first reaches a minimum (corresponding 
to c Si 1.6, assuming an NFW profile) at tacc, when the secondary 
halo first enters the virial radius of the main progenitor, quickly 
followed by a maximum (c Rj 12.0) at tper, when the merging sub- 
halo is at first pericentric approach. A second minimum follows 
next as the accreted material reaches apocenter before relaxing 
to equilibrium. Recent accretion events can clearly bias concen- 
tration estimates unless steps are taken to ensure that halos are 
relaxed. 
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Figure 3. Bottom panels: Halo concentration, as measured by the virial-to-half mass radius ratio, r2oo/'"h; versus half-mass formation 
redshift, z^, for halos in the mass range 3.49 < log M2oo/lO^''^~^-'^0 < 3.79 and resolved with at least 5000 particles in all Millennium 
Simulations (see shaded area in Fig. [ij . Tickmarks on the right indicate concentration values derived assuming that halos follow NFW 
profiles. The different panels show halos identified at three different redshifts: z-^^ = 0, 1, and 3. Grey dots indicate all halos; orange dots 
correspond to halos that pass the relaxation criteria specified in Sec. 12.21 The median concentration of the full halo sample is shown by 
the blue connected circles; thin lines delineate the 25th and 75th percentiles. The shaded area indicate halos whose formation lookback 
time since identification, t^ = tib(^h) ~ilb(^id)i is between one quarter and three-quarter crossing times, defined as tcross = 2r20o/V20o. 
Halos "formed" more than one crossing time ago (to the right of the dot-dashed vertical line) are in general close to virial equilibrium; 
for those halos earlier collapse translates into higher concentration. The dependence of concentration on formation redshift becomes 
non-monotonic when halos that formed less than one crossing time ago are considered, signalling different out-of-equilibrium stages in 
the virialization process. See text for discussion. Top panels: The distribution of formation redshifts for all halos shown in the bottom 
panels (in blue) and for those satisfying the relaxation criteria (orange). The fraction of relaxed halos in this mass range, f^^i, decreases 
with increasing redshift. Very few halos formed less than one crossing time ago pass the relaxation criteria. 



acteristic clustering mass (M.) at that epoch, which is usu- 
ally defined by the condition v{M,) = 1. Tickmarks along 
the top of Fig. [1] list the values of v corresponding to the 
virial masses listed along the bottom axes. At z = even the 
MS-XXL can only probe the regime v <?,, whereas at z = 3 
massive halos trace peaks as rare aav^ 5. These objects are 
considerably rarer than those probed bv lPrada et al.l ()201lh 
(z^ 4), reflecting the extremely large volume of MS-XXL. 

Interestingly, the upturn disappears when only "re- 
laxed" halos are considered, a clear indication that out-of- 
equilibrium effects are responsible. This seems at odds with 
iPrada et al.l (|201ll ). who argue that the upturn remains even 
when unrelaxed halos are excluded; however, they use less 
restrictive relaxation criteria than the ones we adopt here. 
In particular, Prada et al use doff < 0.1 and 2T/\U\ < 1.5, 
compared with our adopted values of 0. 07 and 1.3, respec - 
tively, which are based on the work of iNeto et al.l (|2007l ). 
These relatively small differences in the criteria have a large 
effect on the massive halo sample. For example, at 2 = 3 
fewer than ~ 3% of MS-XXL halos in the sample satisfy our 
offset and virial ratio criteria; on the other hand, 61% of 
them pass the less restrictive criteria adopted by Prada et 
al. 



If the upturn is driven by unrelaxed halos then we 
should be able to get some clues to its origin by considering 
how concentration estimates vary as halos virialize. We turn 
our attention to this point next. 



3.2 Concentration estimates and virialization 

Massive halos formed hierarchically assemble late and they 
are therefore expected to have accreted a large fraction of 
their mass in a short period of time, mainly in the form 
of major mergers. As a merger proceeds, oscillations in the 
mass profile drive fluctuations in the concentration as the 
system relaxes to equilibrium. 

We illustrate this in Fig. [21 where we show the evolution 
of a dark matter halo that undergoes a late-time merger. 
The system chosen fo r this illustration is halo Aq-F from 
the Aquarius Project (|Springel et al.|[2008l V The top panel 
of Fig. [2] shows the evolution of the virial mass of the most 
massive progenitor, normalized to the mass at 2 = 0, as well 
as that of r2oo/''h. The bottom two panels monitor the dy- 
namical relaxation criteria introduced above; when at least 
one of them fails (i.e., strides outside the shaded areas) the 
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Figure 4. Fraction of halos identified at three different redshifts in different stages of relaxation shown as a function of virial mass. The top 
curves correspond to halos that "formed" less than one crossing time ago; i.e., those where t-^ = tib(2h) ~ilb(zid) < tcross = 2r2oo/V2oo. 
These are systems too dynamically young to bo in virial equilibrium. The bottom curves in each panel show systems whose formation 
times suggest that they should have highcr-than-averagc concentrations; i.e., those with 0.25 < th/tcross < 0.75 (see shaded area in 
Fig.O. Different colours are used to indicate results from MS-I, MS-II, and MS-XXL, as indicated in the legend. 



halo would be deemed unrelaxed. We use larger symbols to 
indicate the times along the evolution when this happens. 

Fig. m shows that Aq-F undergoes a major merger with 
a massive secondary halo ("subhalo", for short) at z ~ 0.4. 
Different evolutionary stages are easily identified: the merger 
begins at Zacc = 0.43, when the subhalo first enters the virial 
radius of Aq-F; it reaches pericenter shortly thereafter at 
Zpcr = 0.35. The subhalo is largely disrupted then, and the 
remnant halo bounces once more before quickly virializing 
afterward. 

The effect of these oscillations on the concentration 
are shown by the solid squares that track the evolution of 
^■200 /^'h in the top panel of Fig. [21 Concentrations first drop 
at tacc, reach a maximum at tpcr, and decrease again as the 
subhalo reaches its second apocenter before settling down. 
Since first accretion, it takes at least one crossing time (de- 
fined as fcross = 2r2oo/V^20o) for the system to virialize. The 
arrow labelled t^d indicates the moment when one crossing 
time has elapsed since accretion. 

For a system like Aq-F that nearly doubles its virial 
mass in a merger, the oscillations have large amplitudes: 
roughly a factor of ~ 2 in r2oo/''h, which translates into a 
fluctuation of almost one order of magnitude in concentra- 
tion. Fig. [2] also shows that the dynamical relaxation cri- 
teria introduced above can be effective at identifying sys- 
tems undergoing these transitions. Unless such criteria are 
introduced, the average concentration of a halo sample that 
includes a large fraction of recently-assembled systems may 
differ significantly from that of their equilibrium counter- 
parts. 



3.3 Concentration and formation redshift 

Guided by the results discussed in the previous subsection, 
we explore in Fig. [3] the dependence of concentration on 
formation redshift. Because halos of different mass collapse 
on average at different times, we choose halos in the narrow 
mass range, 3.49 < log(M2oo/10^°/i"^Mo) < 3.79, indicated 



by the shaded vertical band in Fig. [T] This is close to the 
characteristic clustering mass at 2: = (1.35 < u < 1.55), 
but at 2 = 3 corresponds to the rarest, most massive systems 
present in the MS-XXL simulation (4.12 <u < 4.73). 

The three panels of Fig. [3] show results for halos identi- 
fied at 2id = 0, 1, and 3. The bottom panels show the half- 
mass radius concentration estimate, r2oo/rh., as a function 
of the half-mass formation redshift, Zh_, defined as the time 
when the most massive progenitor first reaches one half of 
the virial mass of the halo at zij. Tickmarks on the right in- 
dicate the corresponding NFW-based concentrations. Con- 
nected symbols trace the median concentration; thin lines 
indicate the 25th and 75th percentiles. 

The bottom panels of Fig.[3]show that, as expected, ha- 
los that form earlier tend to be more concentrated. However, 
this holds only for halos that are "old" enough at the time of 
identification to have had a chance to relax to equilibrium: 
the concentration of halos formed less than one crossing- 
time ago (i.e., = fib(2h) - tib(2id) < icross) depends non- 
monotonically on Zh- (We use tib to denote cosmic lookback 
times since z = 0.) The median r2oo/fh trend reverts and 
has a pronounced maximum when ~ fcross/2. The same 
behaviour is observed at all three redshifts, and reflects the 
oscillations in the structure of halos that have recently ac- 
creted a large fraction of their mass as they relax to equi- 
librium. Halos formed in the period 0.25 < th Across < 0.75 
(highlighted by the shaded area in Fig. [3]) are caught at Zid 
during a transient state of maximal concentration. 

In general, the more massive the halo the larger the frac- 
tion of systems that have assembled recently and that are 
therefore out of equilibrium. This is shown in the top panels 
of Fig. [3] where the histograms indicate the distribution of 
formation redshifts of all clusters (in blue) and of relaxed 
clusters (in orange). Clearly, very few halos that have as- 
sembled less than one crossing time ago pass the relaxation 
criteria. Indeed, at z = 3, when halos in this mass range are 
very rare the majority are less than one crossing time "old" 
and fewer than 3% of them are deemed to be in equilibrium. 
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In particular, note that at z = 3 a large fraction of systems 
are at the maximally-compressed stage that corresponds to 
the first pericentric passage of the accreted material, thus 
increasing the average concentration. 

This suggests that the aforementioned upturn in con- 
centration at high masses is due to the increasing preva- 
lence of systems caught during this particular evolutionary 
stage. We show this explicitly in Fig. |4l where we plot, as 
a function of virial mass, the fraction of halos with look- 
back formation times satisfying th < tcross (solid lines) and 
0.25 < ih/tcross < 0.75 (dashed lines). The former character- 
izes halos formed too recently to have reached equilibrium; 
the latter indicates the fraction likely found in a temporary, 
highly-concentrated state. Although recently-assembled sys- 
tems of any mass are susceptible to this effect, the increasing 
fraction of systems in the latter stage with halo mass (or 
provides clear support for our interpretation of the upturn 
as driven by systems passing through a short-lived, highly- 
concentrated stage of their evolution. 

Finally, Fig. |4] makes clear that the fraction of clusters 
out of equilibrium is not just a function of v, but also of red- 
shift. For example, at Zid = 0, about 12% of ^ ~ 2.5 systems 
were formed less than one crossing time ago. At fixed v, this 
fraction increases to 42% and 69% at Zid = 1 and = 3, 
respectively. This increase in the fraction of " unrel axed" 
systems at given i/ explains why IPrada et al.l (|201lh find 
that the median concentration, at given i/, increases with 
redshift (see their Fig. 7). The origin of this effect may be 
traced to the mass-dependent shape of the power spectrum. 
At fixed u, the higher the redshift the lower the halo mass 
and the shallower the slope of the mass fiuctuation spec- 
trum, d log a/ d log M. Shallower spe ctra imply faster struc- 
ture growth (jEfstathiou et allll988l ). implying that, for a 
given V, clusters at high-redshift are being assembled com- 
parably "faster" than their lower-redshift counterparts. This 
increases the out-of-equilibrium fraction with redshift and 
pushes the median concentration even higher. 



4 SUMMARY AND CONCLUSIONS 

We have analyzed dark matter halos identified at different 
redshifts in the Millennium Simulation series (MS-I, MS- 
II, and MS-XXL) to investigate the origin of the recently- 
reported upturn in the mass-concentration relation of rare, 
massive halos. We estimate concentrations from the ratio 
of virial-to-half mass radius, r2oo/'"h, assuming that halos 
follow an NFW profile. Our analysis confirms that, as a 
function of the dimensionless mass parameter u{M,z) = 
ScTit{z) /(j{M, z), the median concentration declines gradu- 
ally until it reaches a minimum at ~ 3, then increases to- 
ward higher values of i'. This upturn, however, is not present 
when only dynamically-relaxed halos are retained for anal- 
ysis, a clear indication that out-of-equilibrium effects drive 
the non-monotonic behaviour of the mass-concentration re- 
lation. 

Further inspection demonstrates that the upturn is due 
to systems caught in a state of "maximal compression" coin- 
cident with the first pericentric passage of recently-accreted 
material. This affects primarily halos that have just accreted 
a large fraction of their mass, and especially those where the 
lookback time to formation coincides approximately with 



Table 1. Halo sample considered in our analysis. Njjaios lists the 
total number of halos with N200 ^ 5000 in each simulation; Njoi 
the total number of "relaxed" halos (see Section 12.211 : and f^^i is 
the relaxed halo fraction. 



Run 




Nhalo 




/rcl 


MSII 





55,152 


42,395 


0.803 




1 


64,650 


35,359 


0.547 




3 


48,321 


10,412 


0.215 


MSI 





83,196 


57,782 


0.695 




1 


64,907 


22,560 


0.348 




3 


6,560 


402 


0.061 


XXL 





2,101,739 


598,603 


0.285 




1 


910,043 


102,475 


0.113 




3 


4,799 


112 


0.023 



the time it takes to bring material from the virial radius of 
the halo to the center, roughly r2oo/V20o- The fraction of 
systems in this transient stage increases systematically with 
halo mass, causing the upturn in the M{c) relation that be- 
comes manifest at extremely high masses, where maximally- 
compressed halos dominate. 

The non-monotonic nature of the mass-concentration 
relation illustrates the fact that, due to their recent assem- 
bly, the vast majority of very massive galaxy clusters are 
generally expected to be out of dynamical equilibrium. It fur- 
ther warns about the applicability of simple power-law scal- 
ings for various observational mass proxies, such as velocity 
dispersion. X-ray luminosity/temperature, or SZ decrement, 
especially when applied to the rarest, most massive clusters. 
Any sample of such systems is likely to be dominated by sys- 
tems currently out-of equilibrium, implying that deviations 
from simple power-laws are to be generally expected. These 
deviations lead to substantial, correlated scatter between 
different mass estimators that should scale in non-trivial 
ways with cluster mass and redshift. Our findings suggest 
that the dynamical state of galaxy clusters must be care- 
fully and explicitly taken into account before cluster studies 
may be used as reliable tools of precision cosmology. 

This paper has been typeset from a T^jiX/ IM]eX file prepared 
by the author. 
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